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ABSTRACT. The opening rate of voltage- 
gated potassium ion channels exhibits a char- 
acteristic, knee-like turnover where the com- 
mon exponential voltage-dependence changes 
suddenly into a linear one. An explanation of 
this puzzling crossover is put forward in terms 
of a stochastic first passage time analysis. The 
theory predicts that the exponential voltage- 
dependence correlates with the exponential dis- 
tribution of closed residence times. This fea- 
ture occurs at large negative voltages when the 
channel is predominantly closed. In contrast, 
the linear part of voltage-dependence emerges 
together with a non-exponential distribution of 
closed dwelling times with increasing voltage, 
yielding a large opening rate. Depending on 
the parameter set, the closed-time distribution 
displays a power law behavior which extends 
over several decades. 



Introduction 

Voltage-dependent ion channels of biological mem- 
branes are formed by pore-like single proteins which 
poke through the cell membrane. They provide the 
conducting pathways for the ions of specific sorts (Q,|2|) . 
Such potassium K + and sodium Na + channels partic- 
ipate in many important processes occurring in liv- 
ing cells. For example, these are crucial for the phe- 
nomenon of neural excitability (^). 

Two features are important for the biological func- 
tion of these naturally occurring nanotubes. First, 
they either are dwelling in open conformations, al- 
lowing for the ion flow to pass through, or are rest- 
ing in closed, non-conducting conformations. Between 
these two conformation types the ion channel under- 
goes spontaneous, temperature driven transitions - the 
so-called gating dynamics - which can be character- 
ized by the residence time distributions of open, f (t), 
and closed, f c (t), states, respectively. The mean open 
and closed residence times, (T ( c )) := f Q tf < c \(t)dt 
are prominent quantifiers of the gating dynamics. In 



particular, they determine the mean opening (closing) 
rates fc„( c ) := (T c r -\) . The second important feature 
refers to the fact that the gating dynamics is voltage- 
dependent. This provides a mechanism for a mutual 
coupling among otherwise independent ion channels. 
This very mechanism is realized through the common 
membrane potential. Both ingredients are central for 
the seminal model of neuronal activity put forward by 
Hodgkin and Huxley in 1952 

The dichotomous character of gating transitions 
yields a bistable dynamics of the Kramers type (f|). 
Therefore, a priori one expects that both, the opening 
and the closing gating rates will expose an exponen- 
tial, Arrhenius-like dependence on voltage and temper- 
ature. Indeed, the closing rate of many K + channels 
follows such a pattern; in clear contrast, however, the 
opening rate usually does not. To explain the exper- 
imental voltage-dependence of the relaxation time of 
the potassium current for a giant squid axon Hodgkin 
and Huxley (|J) postulated that the gating behavior of 
a potassium channel is determined by four indepen- 
dent voltage-sensitive gates, each of which undergoes 
a two-state Markov dynamics with a form (|^,|^) 

a c (V-V c ) 



k (V) 



[1] 



1 - exp[-6 c (^ - Vo)] 

for the opening rate, which is commonly used in neuro- 
physiology. In Eq. |l|, a c ,b c ,V c are some experimental 
parameters. Notwithstanding that in their work (|^) 
this kind of dependence has been used for a single gate, 
the opening rate of the whole K + channel can also be 
fitted by Eq. [j], see e.g. in (||). The same modeling for 
a whole channel is used also for dendritic K + channels 
in neocortical pyramidal neurons @ . 

Note that in Eq. [I] the voltage-dependence of the 
opening rate changes in a knee-like manner from an 
exponential behavior into a linear one, cf. Fig. 1. This 
typical, experimentally observed behavior of delayed 
rectifier K + channels presently lacks an explanation 
in physical terms. A qualitative explanation of this 
gating dynamics has briefly been mentioned in recent 
work (||). However, a definite analysis leading to the 
functional form in Eq. [j] is not available. A first main 
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Figure 1: Dependence of opening (k ) and closing 
(k c ) rates on voltage for a Shaker IR K + channel 
from (0) at T = 18°C. The opening rate is described 
by Eq. (1) with the following parameters @: a c = 
0.03 msec/mV, b c = 0.8 mV" 1 and V c = -46 mV. 
The closing rate is given by k c — 0.015 exp(— 0.038F) 
msec -1 (V in mV) Inset shows the same depen- 

dencies in the semi-logarithmic scale. 



objective of the present work is to fill this gap, and, 
moreover, to provide additional insight into the voltage 
behavior of Eq. |l] within an exactly solvable stochastic 
Fokker-Planck-Kramers model. 

The ion current recordings made on the level of sin- 
gle ion channels (||) reveal yet another unresolved, in- 
teresting aspect of the gating dynamics. Namely, the 
distribution of closed residence times of many channels 
is not exponential. In particular, it has been shown by 
Liebovitch et al. (||) that the closed residence time dis- 
tribution f c {t) in a rabbit corneal endothelium channel 
can be reasonably fitted by a stretched exponential 
with only two parameters. This result initiated the 
construction of the so-called fractal model of ion chan- 
nel gating Other channels, e.g., K + channels 
in neuroblastoma x glioma (NG 108-15) cells exhibit 
a power-law scaling behavior as well, i.e. f c (t) oc t~ a 
with a = 3/2 (11). To explain this type of fractal- 
like behavior Millhauser et al. ( |l2|) proposed a one- 
dimensional diffusion model. Similar power laws with 
a > 3/2 have also been reported and several vari- 
ations of diffusion theory have been introduced to ex- 
plain the gating behavior of different channels (fl4| [l6|) . 

The observed non-exponential behavior can be fitted 
by a finite sum of exponentials; consequently, it can 
alternatively be explained with a corresponding dis- 
crete Markovian scheme (f[l]). These discrete Marko- 
vian models have proven their usefulness in many cases 
(|l7|). Nevertheless, such an approach presents a fitting 
procedure; as such it is intimately connected with the 



Figure 2: Gating dynamics as an activated diffusion 
on a complex free energy landscape. Two global 
minima correspond to open and closed macrocon- 
formations. One assumes a large number of quasi- 
degenerated (within ksT) and voltage-independent 
closed substates separated from the open conforma- 
tion by a voltage-dependent potential barrier. This 
idea is sketched by a simplified model of the Fokker- 
Planck-Kramers type, and by a discrete model with 
open (O), closed (C) and inactivated (I) states. 



danger of a proliferation of parameters. In particular, 
kinetic schemes containing as many as 14 structurally 
unidentified closed substates have been proposed ([l8]) . 

An important lesson to be learned from the detailed 
studies of a simple protein - myoglobin - by Frauen- 
felder et al. ( [T^ ) is that proteins exist in a huge num- 
ber of quasi-degenerated, microscopic substates, cor- 
responding to a single macroscopic conformation, cf. 
Fig. 2. It is thus conceivable that at room tempera- 
tures the ion channel dwells in a huge number of almost 
degenerated (within kgT) conformational substates. 
Both the fractal and diffusion models of the ion chan- 
nel gating have been inspired by this crucial property 
of proteins. We conjecture that the ultimate theory of 
the ion channel gating must take this property into ac- 
count. This program requires a compromise between 
Markovian discrete state models and a continuous dif- 
fusion model. This can be achieved by a Kramers type 
theory (f|,||). The discrete Markov models can then be 
considered as a limiting case of more general Kramers 
type approach (||). 

Theoretical modeling 

The complex structure of the multi-dimensional con- 
formational space of proteins implies an intricate ki- 
netics despite an apparently simple bistability that is 
observed (JL9|). Two popular theoretical approaches 
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have been developed to cope with this complexity. A 
first one uses a simple bistable dynamics as a basis. 
To model the complexity of the observed kinetics this 
dynamics is amended by using an additional stochas- 
tic time-dependence of the energy profile, or kinetic 
constants. Such an approach is nowadays commonly 
known under the label of "fluctuating barriers" (20). 
Alternatively, one can attempt to model the complex- 
ity of the energy profile itself in the simplest possible 
way. Our strategy is to find such a minimal model of 
the second kind which does allow for a rigorous analy- 
sis and does reproduce some nontrivial features of the 
gating dynamics. 

Let us assume that the conformational stochastic 
dynamics between the open and closed states can be 
described in terms of a one-dimensional reaction co- 
ordinate dynamics x(t) in a conformational poten- 
tial U(x), Figs. 2,3. Since the distribution of open 
residence time intervals assumes typically a single- 
exponential (Q), in the following we rather shall fo- 
cus on the behavior of the closed residence time in- 
tervals. In order to evaluate the distribution of closed 
residence time intervals it suffices to restrict our anal- 
ysis to the subspace of closed states by putting an 
absorbing boundary at the interface, x = 6, between 
the closed and open conformations, see Fig. 3. We 
next assume that the gating dynamics is governed by 
two gates: an inactivation gate and an activation gate. 
The inactivation gate corresponds to the manifold of 
voltage-independent closed substates. It is associated 
with the flat part, — L < x < 0, of the potential U(x) 
in Fig. 3. In this respect, our modeling resembles that 
in ( f2lj ) . The mechanism of inactivation in potassium 
channels is quite sophisticated and presently not to- 
tally established (Q) . It is well known that inactivation 
can occur on quite different time scales ([!]). The role of 
a fast inactivation gate in Shaker K + channels is taken 
over by the channel's extended N-terminus which is ca- 
pable to plug the channel's pore from the cytosol part 
while diffusing towards the pore center (22). The slow 
inactivation apparently is due to a conformational nar- 
rowing of the channel pore in the region of selectivity 
filter (@). In both cases, no net gating charge translo- 
cation occurs and the inactivation process does not 
depend on voltage. When the inactivating plug is out- 
side of the pore, or the selectivity filter is open (x > 
in Fig. 3) the channel can open only if the activation 
barrier is overcome. 

The dynamics of the activation gate occurs on the 
linear part of the ramp of the potential U(x); i.e. on 
< x < b in Fig. 3, like in |[|). Note that for 
< x < b, the inactivating plug diffuses outside of the 
channel's pore and the selectivity filter is open. Dur- 
ing the activation step a gating charge q moves across 
the membrane, this feature renders the overall gating 
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Figure 3: Studied model and its diffusion counterpart 

dynamics voltage-dependent. The channel opens when 
the reaction coordinate reaches the location x — b in 
Fig. 3. This fact is accounted for by putting an ab- 
sorbing boundary condition at x — b. Moreover, the 
channel closes immediately when the inactivation gate 
closes (x < 0), or when the activation gate closes. To 
account for this behavior in extracting the closed res- 
idence time distribution we assume that the channel 
is reset into the state x = after each closure (see 
below) . 

The diffusional motion of the inactivated gate is re- 
stricted in conformational space. We characterize this 
fact by the introduction of a conformational diffusion 
length L (Fig. 3) and the diffusion constant D ~ ksT 
that are combined into a single parameter - the con- 
formational diffusion time 



t d = L 2 /D 



[2] 



This quantity constitutes an essential parameter for 
the theory. We assume that the activation barrier 
height Uq is linearly proportional to the voltage bias 
V jl^), i.e. in terms of the gating charge q we have 



U 



-q{V-V c ). 



[3] 



Moreover, Uq is positive for negative voltages, i.e. for 

V < V c , vanishes at V — V c , and becomes negative for 

V > V c . Thus, for V > V c the channel "slips" in its 
open state, rather than overcomes a barrier. In addi- 
tion, the fraction £ of the voltage-dependent substates 
in the whole manifold of the closed states should be 
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very small, implying that 



With the parameter identifications 



e = 6/x«i. 



[4] 



Analytical solution. The corresponding Fokker- 
Planck equation for the probability density of closed 
states P(x,t) reads 

dP(x,t) d ( d n 8U(x)\ . , 

where (3 — \/{kgT). In order to find the distribution 
of closed residence times f c (t), we solve Eq. ||]with the 
initial condition P(cc,0) = 5(x), in combination with 
a reflecting boundary condition dP ^'*' ) \ x =-l = 0, and 
an absorbing boundary condition, P(x, t)\ x =b = (|4|). 
The closed residence time distribution then follows as 



fc(t) 



d<S> c (t) 
dt 



[6] 



where & c (t) = f_ L P(x, t)dx is the survival probability 
of the closed state. 

By use of the standard Laplace transform method 
we arrive at the following exact solution: 



Us) 



A(s) 

B( s y 



where 



[7] 

[8] 
[9] 



A(s) = exp(-/3C/ /2) v //3 2 i7 2 + 4^r flS 

B(s) = y/ /PU* + 4?t d s 

x cosh [\jf3 2 Ul + 4£ 2 t d s/2) 

+ ^2^rDstanh ^/tds — f3Utjj 

x sinh (J (3 2 U$ + 4e 2 r fl s/2) . 

The explicit result in fr|-^ allows one to find all mo- 
ments of the closed residence time distribution. In 
particular, the mean closed residence time (T c ) = 
lim s ^ [l - f c {s)]/s reads 



(T c ) = t d £ 



PU (e 



PU 



/3Ua 



[10] 



This very same result ^ can be obtained alternatively 
if we invoke the well-known relation for the mean first- 
passage time (T c ) = i ^dxe^^ J^dye-W^ (g). 
This alternative scheme provides a successful validity 
check for our analytical solution in [7||^. 

Elucidation of the voltage dependence in Eq. 
[l| Upon observing the condition^ Eq. 10 by use of ||| 
reads in leading order of £ 



1 

W) 



(3q 



V-V R 



fill 



and 



[12] 



[13] 



the result in 
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precisely coincides with Eq. |l|. The 
fact that our novel approach yields the puzzling volt- 
age dependence in Eq. l] constitutes a first prime result 
of this work. 

Let us next estimate the model parameters for a 
Shaker IR K + channel from Ref. (^). In ([]]), the 
voltage-dependence of k (V) at T = 18 °C has been 
parameterized by Eq. [j] with the parameters given in 
the caption of Fig. 1. Then, from Eq. 12 the gating 
charge can be estimated as q ~ 20e (e is the positive 
valued, elementary charge). As to the diffusion time 
t d , we speculate that it corresponds to the time scale 
of inactivation; the latter is in the range of seconds 
and larger (Q). Therefore, we use tjj = 1 sec as a 
lower bound for our estimate. The fraction of voltage- 
dependent states £ is then extracted from Eq. |l^ to 
yield, £ sa 0.0267. This value, indeed, is rather small 
and thus proves our finding in Eq. |ll] to be consistent. 

Analysis for the closed residence time dis- 
tribution. The exact results in Eqs. appear 
rather entangled. To extract the behavior in real time 
one needs to invert the Laplace transform numerically. 
With £ << 1, however, Eqs. are formally reduced 
to 



U*) 



1 



1 + (k To) ^-Jtds tanh */tdS 



[14] 



This prominent leading order result can be inverted 
analytically in terms of an infinite sum of exponentials, 
yielding: 



Ut) = X! Cn ^ n cxp(-A„£), 



[15] 



where the rate constants < Ai < A2 < ... are solu- 
tions of the transcendental equation 



tan yX, 



k T D 

vAnTD 



[16] 



and the expansion coefficients c„, respectively, are 
given by 



[17] 



1 + k Q TD + X n /k 



Note from Eq. |^ that the set c n is normalized to unity, 
The analytical approximation, Eqs. 15-17, is com- 



pared in Fig. 4 with the precise numerical inversion 
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Figure 4: Closed residence time distribution for a 
diffusion-limited case. The exact numerical result (full 
line) is compared with the analytical approximation in 
Eqs. (15)-(17) (broken line). The latter one coincides 
with the exact solution of the diffusion model by Mill- 
hauser et al. in the scaling limit. 

of the exact Laplace transform in Eqs. ^-^|. The nu- 
merical inversion has been performed with the Ste- 
hfest algorithm (|2^). As can be deduced from Fig. 
(4), for t > 10 msec the agreement is very good in- 
deed. A serious discrepancy occurs only in the range 
0.01 msec < t < 0.1 msec which lies outside the range 
of the patch clamp experiments (t > 0.1 msec). More- 
over, the agreement is improving with increasing td 
(not shown). 

Origin of the power law distribution. The fea- 
tures displayed by the closed residence time distri- 
bution f c (t) depend sensitively on the applied volt- 
age V. When V > V c , e.g. V = —45 mV, as 
in Fig. 4, the activation barrier towards the chan- 
nel opening disappears and the opening dynamics be- 
comes diffusion-limited. In this case, the diffusion time 
td = 1 sec largely exceeds the mean closed residence 
time (T c ) « 18.4 msec. Put differently, t d > (T c ) and 
the closed residence time distribution exhibits an in- 
tricate behavior with three distinct regions, see in Fig. 
4. Most importantly, for the intermediate time scale 



(T C ) 2 /T D « t « T D 



[18] 



we find from Eq. 14 (by considering the limit td — > 
00) that the closed residence time distribution obeys a 
power law; reading 



fc(t) 



1 



2(7TT D ) 1 /2/U3/2- 



[19] 



This type of behavior is clearly detectable in Fig. 4 
where it covers about two decades of time. As follows 
from Eq. 18, an increase of td by one order of magni- 
tude (while keeping (T c ) fixed) extends the power law 



region by two orders of magnitude. This conclusion 
is fully confirmed by our numerics (not shown). This 
power law dependence, which extends over four orders 
of magnitude, has been seen experimentally for a K + 
channel in NG 108-15 cells ([Tl"|). On the contrary, for 
channels, where td is smaller, the power law region 18 
shrinks and eventually disappears, whereas the mean 
opening rate defined via Eq. ^ still exhibits a steep 
dependence on the voltage. Thus, our model is capable 
to describe for different channels both, the emergence 
of power law as well as its absence. 

On the time scale t > td the discussed power law 
distribution crosses over into the exponential tail; the 
latter is fully described by the smallest exponent Ai in 
Eq.15, i.e., by 



f c (t) « C1A1 exp[-Aii] 



[20] 



This feature is clearly manifest in Fig. 4. The tran- 
sition towards the exponential tail in the closed res- 
idence time-interval distribution can be used to es- 
timate the diffusion time td on pure experimental 
grounds! 

Finally, let us consider the opposite limit, td <C 
(T c ), for V <C V c . For the considered set of param- 
eters this occurs, e.g., for V = —55 mV when the 
channel is predominantly closed. Then, the diffusion 
step in the opening becomes negligible and in the ex- 
perimentally relevant range of closed residence times, 
defined by (T c ) , the corresponding distribution can be 
approximated by a single exponential 20 . A perturba- 
tion theory in Eq. 16 yields Ai ss fc D (l — (fc T£))/3). For 
the used parameters we have Ai « 0.96fc o and, from 
Eq. [l7|, ci ps 0.95. This result is in a perfect agree- 
ment with the precise numerics obtained from Eqs. 
0-^. Thus, the distribution of closed residence times 
is single-exponential to a very good degree. Conse- 
quently, one and the same channel can exhibit both, 
an exponential and a power-law distribution of closed 
residence times, as a function of the applied trans- 
membrane voltage. With an increase of td the voltage 
range of the exponential behavior shifts towards more 
negative voltages, V < V c , and vice versa. 

Reduction to a diffusion model. Let us relate 
our model to that introduced previously by Millhauser 
et al. (|lj). The latter one is depicted with the lower 
part in Fig. 3. It assumes a discrete number N of 
closed substates with the same energy. The gating 
particle jumps with the equal forward and backward 
rates k between the adjacent states which are occupied 
with probabilities p n (t). At the right edge of the chain 
of closed states the ion channel undergoes transition 
into the open state with the voltage-dependent rate 
constant 7. To calculate the closed residence time dis- 
tribution f c {t) one assumes po(0) = 1, p n ^o(0) = and 
d® c {t)/dt = -7Po(t), where $ c (i) = En=Q iV Pn(*) is 
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the survival probability (]l 2j.|l-'j|) . 

We consider the continuous diffusion variant of this 
model (24) in a scaling limit: wc put Aa; — > 0, k — ► 
00,7 — » 00, A — > 00 keeping the diffusion length i = 
AAx, the diffusion constant Z? = fc(Aa;) 2 , and the 
constant fc G = 7/iV all finite. The latter one has the 
meaning of mean opening rate, see below. Note that in 
clear contrast with our approach, the rate parameter 
k a in the diffusion model is of pure phenomenological 
origin. The problem of finding the closed residence 
time distribution is reduced to solving the diffusion 
equation 



dP(x,t) __ D d 2 P(x,t) 



ot 



dx 2 



[21] 



with the initial condition P(x,0) — 5(x — 0_), the 
reflecting boundary condition dP ( x ^ \ X= _ L = and 
the radiation boundary condition (pq) 



dP(x, t) l 
dx 



x=0 



Lk 



P(0,t). 



[22] 



We emphasize that the radiation boundary condition 
22 is not postulated, but is rather derived from the 
original discrete model in the considered scaling limit. 
Using the Laplace transform method we solved this 
problem exactly and obtained the result in Eq. 14. In 
conclusion, our approximate result in Eqs.14-17 pro- 
vides the exact solution of the diffusion model (jl2, 15) 
in the scaling limit! This exact analytical solution is 
obtained here for the first time. Note, however, that 
this so obtained diffusion model is not able to resolve 
the puzzling voltage dependence in Eq. |l|. 



Synopsis and Conclusions 

With this work we put forward a unifying generaliza- 
tion of the diffusion theory of ion channel gating by 
Millhauser et al. (|l2,15). Our novel theory reproduces 
for the first time the functional form of the puzzling 
voltage-dependence in Eq. |l[ The latter has been pos- 
tulated almost fifty years ago in the pioneering paper 
by Hodgkin and Huxley (j^) and is commonly used in 
the neurophysiology up to now. The proposed model 
of the Fokker-Planck-Kramers type explains the ori- 
gin of steep voltage-dependence in Eq. [I] within a 
clear physical picture which seemingly is consistent 
with both our current understanding of the physics of 
proteins and basic experimental facts. Our study fur- 
thermore reveals the connection between the voltage 
dependence of the opening rate and the intricate be- 
havior for the closed residence time distribution in cor- 
responding voltage regimes. A particularly appealing 
feature of our approach is that our model contains only 



four voltage-independent physical parameters: the dif- 
fusion time td , the fraction of voltage -dependent sub- 
states £, the gating charge q and the threshold voltage 
V c . Several experimental findings could be described 
consistently while still other ones call for an experi- 
mental validation. 

In particular, when (i) the activation barrier is very 
high, i.e., V <C V c , the activation step determines com- 
pletely the opening rate: the distribution of closed 
residence times is nearly exponential, as well as the 
voltage-dependence of the opening rate. The chan- 
nel is then predominantly closed. We remark that 
the opening rate should exhibit an exponential depen- 
dence on temperature as well. This conclusion fol- 
lows from Eqs. [Tl]-|l2| and the fact that in accord 
with our model the parameter a c in Eq. |l| is tem- 
perature independent. Indeed, with the diffusion time 
td being inversely proportional to the temperature, 
i.e. t d ~ IAD ~ l/(k B T), yielding a c ~ l/( TD k B T), 



cf. Eq. 13. In contrast, when (ii) the activation 
barrier vanishes, i.e. the voltage shifts towards the 
positive direction, the closed residence time distribu- 
tion becomes non-exponential. On the intermediate 
time scale given in Eq. |l8| this distribution exhibits 
a power law behavior, f c (t) oc £~ 3 / 2 , which crosses 
over into an exponential one at t > tjj. The emer- 
gence of the exponential tail can be used to determine 
the conformational diffusion time td experimentally. 
When (iii) the activation barrier assumes negative val- 
ues at voltages V > V c , our result for the opening rate 
exhibits a linear dependence on voltage and, conse- 
quently, see Eq. 
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it no longer depends on temper- 
ature. The weak temperature dependence will emerge 
however when we renormalize the diffusion coefficient 
D due to the roughness of random energy landscape 
(cf. Fig. 2). Assuming uncorrelated Gaussian disor- 
der one gets D ~ k B Texp{-(6U 2 )/(k B T) 2 ) 
where (SU 2 ) is the mean-squared height of the barrier 
between substates. Then, k D ~ exp(— (SU 2 ) / (k B T) 2 ) , 
and since y (SU 2 ) ~ k B T this non-Arrhenius depen- 
dence is weak at room temperatures. This result has 
a clear thermodynamic interpretation: when the ac- 
tivation barrier vanishes the closed-to-open transition 
is entropy dominated and thus the opening rate will 
only weakly depend on temperature. In accord with 
our model this type of behavior correlates with a non- 
exponential kinetics. 

The temperature behavior of the opening rate (or, 
equivalently, the mean closed time) presents a true 
benchmark result of our theory. The authors are look- 
ing forward to seeing this feature being tested experi- 
mentally 
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